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ABSTRACT 

We study the relation between the star formation intensity of galaxies and the extinc- 
tion by dust of their emitted light. We employ a detailed statistical analysis of Hubble 
Deep Field North data to show a clear positive correlation between the extinction and 
star formation intensity at all epochs from redshift 0.4 to 6.5. The extinction evidently 
increases with time for a given star formation intensity, consistent with the expected 
increase in the metallicity with time. Our observational results are well fitted at all 
epochs by a double power-law model with a fixed shape that simply shifts with red- 
shift. The correlation between the extinction and the star formation intensity can be 
interpreted by combining two other trends: the correlation between the star forma- 
tion rate and the gas content of galaxies, and the evolution of the dust-to-gas ratio 
in galaxies. If we assume that Kennicutt's observed relation for the former is valid at 
each redshift, then our findings imply an interesting variation in the dust-to-gas ratio 
in galaxies within each epoch and with time, and suggest new ways to investigate the 
cosmic evolution of this quantity. 
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1 INTRODUCTION 

In this Letter we investigate the relation between the ex- 
tinction and the star formation rate (SFR) in galaxies at 
various epochs. Understanding this relation is a key for bet- 
ter understanding the physical conditions in the interstellar 
medium (ISM) of galaxies as well as the environments of 
galaxies and their cosmic evolution. The SFR in galaxies has 
been studied extensively in recent years. The mo st widely 
used s tar formation law, originally conjectured by ISchmidtl 
1 19591 ). correlates the SFR in galaxies and their gas density 
as a p ower law. As summarized and integrated bv lKennicuttl 
1 19981 ) ■ many empirical studies have tested and proven the 
effectiveness of the Schmidt law. Now, extinction by dust 
is related to star formation since dust forms from heavy 
elements ejected during stellar mass loss. Also, the dust-to- 
:as ratio is re lated to the metallicity of the interstellar gas 
Schmidt fc Boile r 1993; Kim & Martin 1996). 

It is thus clear that a correlation should exist between 
the SFR and the extinction. Indeed, previous observations 
indicate such a correlation w ithin galaxies, and fo r galax - 
ies at v a rious r edshifts [e.g., Adelberger fc Steidell l|2000l ) ; 
ICalzettil (|200ll ): ICalzetti et aLriiroOSl )]. The goal of this pa- 
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per is to observe this correlation and determine its form for 
a broad range of redshifts, and to demonstrate that our con- 
clusions are not strongly affected by selection effects. In the 
following section we describe the dataset used for the anal- 
yses in this Letter. In section 3 we test the dataset to verify 
that it is statistically valid and can be used for the analy- 
ses to follow. We present our main results and model fits in 
sections 4 and 5, and give our conclusions in section 6. We 
note that we focus on the correlation of extinction with the 
star formation intensity (SFI) of galaxies (which is the SFR 
per unit area), but consider in § 5 also the correlation with 
the total SFR. 



2 THE DATA 

We used the dataset described in detail in Thompson, 
Waymann fc Sto rrie-Lombardi (2001; hereafter TWS) and 
iThompsonI (|2003l ). It is a combination of observations of the 
Hubble Deep Field North (HD F-N) with the Wide Field 
Planetary Camera 2 [WFPC 2; IWiUiams etHI 11996")] and 
with the Near-I nfrared Camera a nd Multi-Object Spectrom- 
eter [NICMOS: [bickinsonl (|2000l )). The main output o f the 
photometric analysis in TWS and iThompsonI (j2003h was 
the most likely combination of redshift, extinction and spec- 
tral energy distribution (SED) for each of the galaxies that 
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passed the selection criteria. However, in this present work 
we utilized not only the most likely combinations but the en- 
tire database of the photometric analysis. The idea is that 
while the uncertainties are large from each photometric fit to 
an individual galaxy, if we include these uncertainties prop- 
erly then the total dataset of thousands of galaxies will have 
strong statistical power for studying galaxy evolution. 

The dataset includes two three-dimensional matrices for 
each of the 1972 galaxies that passed the selection crite- 
ria and were used for the analysis. The matrix dimensions 
are 101 x 51 x 15, corresponding to 101 possible redshift 
values between and 8, 15 extinction values ranging from 
E{B -V) = 0to 1.0, and 51 possible SED templates. The 
first matrix holds the values for each of the redshift, ex- 
tinction and SED template combinations, determined as de- 
scribed in TWS by computing the goodness-of-fit to the six 
observed photometric fiuxes. Note that in computing the x^, 
TWS added a term to the background noise that is 10% of 
the flux, to represent various systematic effects. As discussed 
further in TWS, this term may overestimate the noise some- 
what and lead to an underestimate of the values, produc- 
ing in our analysis below a slightly broader distribution of 
probable parameter values. The second matrix in the dataset 
holds the UV flux at 1500A rest-frame that corresponds 
to each of the redshift, extinction and template combina- 
tions. We used this matrix to determine the SFR that corre- 
sponds to the different combi nations according to the broad- 
band 1500A relation g iven by Madau. Pozzetti. fc DickinsonI 
l|l998h . modified for a lScalol (frggsh IMF: 



Fisoo = 10.0 X 10^'' SFR [Mgyr"^] erg s"^ Hz ^ . (1) 

In addition to the above photometric dataset, for 131 
:alaxi es spectroscopic redshifts are available from ICohenI 
200l[ l. For these galaxies, we adopted the spectroscopic red- 
shift and used the above matrices to determine the likelihood 
of various extinctions and SED templates given the known 
redshift. 



3 TESTING THE DATASET USING 
REDSHIFT ERROR ANALYSIS 

In general, a x^ value is related to the corresponding prob- 
ability through the relation: 



P oc exp 



(-1 



X 



(2) 



We therefore used this relation to convert, for each galaxy, 
the x^ values that correspond to the different combinations 
of redshift, extinction and template, into their correspond- 
ing probabilities. The combinations for each galaxy were 
then divided into two groups, one with redshift values higher 
than the redshift value of the most probable combination, 
and the other with redshift values that are lower than the 
redshift value of the most probable combination. We then 
found 1-a redshift confidence limits upwards and downwards 
by identifying the redshift values at 0.683 of the cumulative 
probability of each of the two groups. We emphasize that a 
full probability distribution of possible parameters is deter- 
mined by fitting to each observed galaxy. Figure [T] compares 
the photometric redshift values with errors to the spectro- 
scopic values for the 131 galaxies for which spectroscopic 
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Figure 1. Photometric redshift versus spectroscopic redshift for 
each galaxy (in a log plot of 1 + ^). The confidence limits of l-tr 
up and 1-cr down for the photometric redshifts were determined 
using the probability distribution of the photomet ric redshift o f 
each galaxy. The spectroscopic data are taken from lCohenI (l200lh 



data are available from ICohenI ||200 



to the plot in figure 3 in iThompso: 



. This plot is similar 
(|2003l ) with the sig- 
nificant difference that here the 1-a confidence limits take 
into account the entire distribution of possible photometric 
redshifts, a distribution that in some cases is very wide or 
multiply-peaked. 

The photometric redshifts and their calculated errors, 
together with the spectroscopic redshift data, can be used 
to test the validity of the photometric analysis dataset. This 
test is important as the analysis in the following sections 
strongly relies on the statistical properties of the x'^ 
trices. The photometric redshifts deviate from the spectro- 
scopic ones, but we now consider how these deviations are 
distributed. If the statistical analysis is valid, then the dis- 
tribution of the deviations between the photometric and 
the spectroscopic redshift should be roughly Gaussian, with 
standard deviation as predicted by the statistical analysis. 
Rather than restricting to a symmetric distribution, we ap- 
ply the separately-determined upper and lower confidence 
limits. Figure [2] shows the normalized relative redshift error 
distribution (upper panel). We define 5^ as the fractional 
redshift error expressed in standard deviations: 



^phot I 



(3) 



where a± is the upper or lower confidence limit which is cho- 
sen according to whether the spectroscopic redshift is higher 
or lower than the photometric one, respectively. In the Fig- 
ure we also show (lower panel) the cumulative distribution 
of all the points, both those with positive and those with 
negative errors. We find that the empirical distribution is 
indeed close to a normal distribution with a standard de- 
viation that agrees well with that predicted by our fitting 
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Figure 2. Upper panel: The normalized distribution of the frac- 
tional redshift error expressed in standard deviations, , defined 
in equation Lower panel: The cumulative distribution of Sz- 
The solid curves in the two panels represent the theoretical Gaus- 
sian and error function curves, respectively. 



method; in particular, the cumulative 1-cr, 2-a and 3-cr val- 
ues are 0.702, 0.893 and 0.969, compared to the theoretical 
expectation of 0.683, 0.954, and 0.997, respectively. 



4 CORRELATING THE EXTINCTION WITH 
THE STAR FORMATION INTENSITY 

The SFI (denoted x) is defined as the SFR per area, ex- 
pressed in solar mass per year per proper kiloparsec squared. 
As described in section 2, various combinations of redshift, 
extinction and template may be associated with each galaxy 
and for each of these combinations the SFI can be deter- 
mined. First the SFR of each combination is determined 
from equation ([1} , and then the area is calculated (for each 
possible redshift) from the number of pixels that the galaxy 
covers in the HDF-N image. The average SFI of the galaxy 
is then the SFR divided by the proper area. We emphasize 
that in the analysis below we use the average extinction and 
SFI of each galaxy, and do not consider individual pixels 
separately. In this work we assume in all cosmological con- 
versions a flat cosmology with h — 0.72, flm = 0.27 and 
= 0.73. 

We now consider our main objective in this Letter which 
is to correlate the extinction and the SFI. As noted in sec- 
tion 3, the probabilities of the different combinations of ex- 
tinction, redshift and template for each galaxy can be calcu- 
lated using the values associated with the combinations 
together with equation ([2]). Having obtained the probabili- 
ties, for each galaxy we statistically sample this distribution 
200 times in a Monte Carlo approach. The data that are used 
in the following analysis are the redshift, extinction and the 
calculated SFI values, for each of the 200 sample combina- 



Figure 3. Upper panel: Plot of the extinction E(B — V) versus 
X for the redshift bin z = 0.7-1.0. The detection threshold is 
apparent as a linear cutoff at the upper-left end of the distribution 
of data points (i.e., at low x or high E{B—V) values). Lower panel: 
Histogram of the x values in the upper panel that correspond to 
E{B — V) = 0.2. Note that the peak lies far from the detection 
limit. 



tions, for each of the 1927 galaxies used in this Letter. This 
Monte-Carlo approach includes statistically all the uncer- 
tainties and degeneracies from the photometric fits of the 
individual galaxies. Figure [3] (upper panel) plots the extinc- 
tion versus x for combinations with redshifts z = 0.7-1.0. 
The binning of the extinction values is apparent. There are 
12,072 points in the plot and their distribution clearly points 
to a correlation between the extinction and the SFI. 

We note that some caution is needed here, since a corre- 
lation between the extinction and the SFI would have been 
introduced into our analysis simply as a selection effect; thus 
we must carefully demonstrate that the observed correlation 
is instead a real effect. To identify an object as a galaxy, the 
selection procedure (as described in TWS) requires a mini- 
mum flux per pixel in the various filter wavebands. At any 
given redshift, objects with very low x values or with high 
extinction values may exist and not be detected. This se- 
lection effect is evident in the upper panel of Figure [3] as 
a lack of galaxies at the low x and high E region (i.e., the 
top-left corner). More specifically, the left boundary of the 
collection of data points is quite linear (as is the case in 
other redshift bins as well) . Since the detectability of galax- 
ies is roughly set by a minimum surface brightness, and the 
observed surface brightness is proportional to the intrinsic 
X of the galaxy times an extinction correction of the form 
exp[—aE{B — V)] for some constant a, then galaxies with 
observed fiux at the detection limit should lie on a straight 
line in E versus \og^Q{x). 

In general, a correlation between the extinction and the 
SFI could be inferred mistakenly due to this selection effect, 
but in Figure |3] this is not the case. The selection effect is 
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relevant only at low x and high E, but figure (3a) also shows 
a clear lack of points at high x and low E, where any galaxies 
would be easily detected, and selection effects are absent. 

To demonstrate this more clearly, the data points at 
each extinction value were tested to verify that most of the 
X values lie far above the detection limit, and thus the mean 
X value (which we use below) is insensitive to selection ef- 
fects. The lower panel of Figure |3] shows a histogram of the 
logj^g(x) values that correspond to E{B — V) = 0.2 in the 
upper panel of the same figure. As can be seen, the detection 
limit in this case is around log]^o(^) ~ "Sj whereas most of 
the data points lie at logj^g(a;) — —1.6 or higher, and the av- 
erage is at -1.1. This behaviour, which we have verified for 
other extinction values and other redshift bins, clearly indi- 
cates that the selection effect is not the main cause of the 
correlation we find between the extinction and SFI. We cau- 
tion though that at the highest redshifts we consider, there 
are fewer observed galaxies and fewer data points particu- 
larly at high extinction, so the detection threshold affects 
the results significantly in this regime. Additional data at 
fainter fluxes would be very useful for checking the reliabil- 
ity of the high-redshift, high-extinction points. 

In order to characterize the correlation between the ex- 
tinction and the SFI as a function of epoch, we divided the 
data into ten redshift bins: 0.4-0.7, 0.7-1.0, 1.0-1.5, 1.5- 
2.0, 2.0-2.5, 2.5-3.0, 3.0-3.5, 3.5-4.5, 4.5-5.5 and 5.5-6.5. 
In each redshift bin, the mean and standard deviation of all 
the logjQ(a;) values that correspond to each speciflc extinc- 
tion value were calculated. We calculated symmetric errors 
here in order to use them below in a minimization. The 
resulting extinction and mean logjQ(a;) values, with 1-a error 
bars in logjg(a;), are plotted in Figure |4] for all redshift bins. 
The correlation between the extinction and SFI is clearly 
evident at all epochs. Note that the errors in x [which is 
proportional to (1 -|- z)*] are smaller in the higher-redshift 
bins since the photometric redshifts are more accurate due 
to the Lyman break (see also figure 1). 

Figure [5] (upper panel) plots all the redshift bins to- 
gether so that the evolution with redshift can be observed 
more easily. It is evident that for a given SFI, the extinction 
decreases with z. This behaviour is consistent with the in - 
crease in the metallicity of the ISM with time l|lnouell2003l) . 
The lower panel of Figure [5] is similar and is based on the 
same analysis as the upper panel except that the extinction 
is plotted against the total SFR and not against the SFI. 
Although the results in the lower panel depend differently 
on the size distribution of galaxies, the same behaviour is 
seen as in the upper panel. We further consider these two 
different correlations in the following section. 



5 FITTING THE CORRELATION 

In the previous section, the relation between the extinction 
and SFI was described empirically. As noted in the introduc- 
tion, this relation is expected from the Schmidt law together 
with the correlation between the gas and dust in the ISM of 

gala xies. 

ISchmidtl (|l959h put forth the hypothesis that the rate 
of star formation in a given region varies as a power of the 
gas density within that region. Observations yield direct esti- 
mates of surface densities and thus, the relation between the 
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Figure 4. Plots of the extinction E{B — V) versus the SFI at 
various redshift bins between z = 0.4 and 6.5. The correlation 
between the extinction and the SFI is clearly observed at all 
these epochs. We show the double power-law fit to the data (solid 
curves). We also compare to a power-law slope of 0.71 (dashed 
curve, z = 0.7-1.0 panel). See section 5 for a discussion of the 
model fits 
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Figure 5. Plot of E{B — V) versus x (upper panel) or versus 
the SFR (lower panel) for all redshift bins. The bins are z = 0.4- 
0.7 (x), 0.7-1.0 (•), 1.0-1.5 (*), 1.5-2.0 (open squares), 2.0-2.5 
(+), 2.5—3.0 (solid triangles), 3.0—3.5 (o), 3.5—4.5 (solid squares), 
4.5—5.5 (o) and 5.5—6.5 (open triangles). 
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SFI X and the gas surface density has been measured over 
a broad range of physical conditions in relatively nearb y 
galaxies, yielding the relation x oc Eg^t (jKennicuttI 1 1998l l . 
The relation between the dust and gas content in the ISM is 
less clear. If a constant dust-to-gas ratio is assumed, then if 
the Kennicutt relation is valid at a given redshift it should 
yield a linear relation between x and E, with a slope of 
1/1.4 = 0.71 in a log-log scale. However, this model is in- 
consistent with the observational results shown in Figure 
To illustrate this, we show in the top-right panel {z — 0.7- 
1) of this figure a line with this slope (dashed curve). It is 
apparent in all the redshift bins that the dependence of E 
on X in the actual data is steeper than this power at low x 
and shallower at high x. Equivalently, the dust-to-gas ratio 
at both ends is lower than at the central x values, assuming 
the lKennicutd (| 19981 ) relation is valid at each redshift. 

The behaviour in Figure |4] suggests a double power- law 
model; the precise form we choose is motivated by the con- 
ventional model for the luminosity function of quasars [e.g., 
IPetersonI ()l997l )]. In addition. Figure [5] suggests a simple 
evolution with redshift as a power of {1 + z), and so we fit 
a single model to all epochs simultaneously: 



E{B -V) = 



Eo{l 



+ 



xs(l + z)"'2 



(4) 



This model consists of six free parameters: Eo - a. normaliza- 
tion parameter, xs - the knee point that separates the two 
slopes, 7i and 72 - the powers of (1 -I- z) that describe the 
redshift evolution of Eo and xs, respectively, and a and /3, 
which represent the asymptotic power-law slopes at low and 
high a;-values, respectively. Note that we excluded E = 
data points from the fit and from Figures |4] and (5] 

We fit the model to all redshift bins by minimizing a 
combined built according to the errors shown in Figure [l] 
The figure shows that the model yields a rather good fit to 
the data at all epochs. The redshift value of the fit shown in 
each panel is the value at the center of the redshift bin. The 
fitted parameters with their 1-a errors (which correspond to 
a unit increase in the value from its minimum value) are: 



-0.65t;^?, log,o(a:s) 



-2.79^ 



71 



0^ 



logio(-Bo) 

72 = 4.14t;;t^, a = 2.03+:42, and /3 = 0.41+;^^. 

The best-fit value 71 = implies a relation of constant 
shape that simply shifts horizontally with redshift. The 
value of the best fit is 111.6 while the number of degrees 
of freedom is 134, which suggests that the errors may be 
somewhat overestimated, which is perhaps related to our 
assuming uncorrelated errors for the various points when 
constructing the for this fit. Figure [5] reveals another 
possible characterization of the correlation by relating the 
extinction to the SFR instead to the SFI. The double power 
law was fitted to these data as well, using the same anal- 
ysis as described above for the case of the SFI. The fitted 
parameter results in this case, with a best-fit value of 



107.1, are: logi(,(£;o) 



71 



-0.39^:1;?, logio(a;s) 



: 2.161 



^, a = 2.17: 



-0.62 



and P = 0.22: 



+ .14 
-.13' 
05 
06- 



The two fits presented here depend on the sizes of galax- 
ies in different ways. While the SFI is a more local quantity, 
the total SFR depends also on the global size of the galaxy. A 
galaxy's overall size and mass affects star formation through 
global feedback effects. In particular, supernova feedback is 
efficient only in relatively small galaxies which results in var- 



ious observe d correlations in the proper ties of galaxies at low 
redshift [e.g.. lDekel fc BirnboimI l|2005h ]. Additional data for 
a larger number of galaxies is required to fully explore the 
effect of the size of a galaxy on its average extinction value. 



6 CONCLUSIONS 

We have shown that a clear correlation exists between the 
extinction of galaxies and their SFI. This correlation is ap- 
parent at various epochs ranging from z—0.4 to z=6.5 and 
is in general agreement with an increasing metallicity in the 
ISM with time. We have also shown that a similar correla- 
tion remains when the galaxies' SFR is considered instead of 
their SFI, a case in which global effects such as supernovae 
feedback may affect the behaviour. We have fitted a double 
power-law model to the results and shown that it fits the 
data well. We have suggested that this correlation is a natu- 
ral consequence of the Schmidt law, which relates the SFI of 
galaxies with their gas content in the ISM, together with a 
relation that exists between the gas and dust content in the 
ISM. The exact shape of the correlation should be further 
investigated both to further eliminate selection effects and 
to shed more light on the cosmic evolution of gas and dust 
in galaxies. 
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